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On-chip magnets can be used to implement relatively large local magnetic field gradients in na- 
noelectronic circuits. Such field gradients provide possibilities for all-electrical control of electron 
spin-qubits where important coupling constants depend crucially on the detailed field distribution. 

We present a double quantum dot (QD) hybrid device laterally defined in a GaAs / AlGaAs het¬ 
erostructure which incorporates two single domain nanomagnets. They have appreciably different 
coercive fields which allows us to realize four distinct configurations of the local inhomogeneous 
field distribution. We perform dc transport spectroscopy in the Pauli-spin blockade regime as well 
as electric-dipole-induced spin resonance (EDSR) measurements to explore our hybrid nanodevice. 
Characterizing the two nanomagnets we find excellent agreement with numerical simulations. By 
comparing the EDSR measurements with a second double QD incorporating just one nanomagnet 
we reveal an important advantage of having one magnet per QD: It facilitates strong held gradients 
in each QD and allows to control the electron spins individually for instance in an EDSR experi¬ 
ment. With just one single domain nanomagnet and common QD geometries EDSR can likely be 
performed only in one QD. 

PACS numbers: 73.63.-b, 03.67.-a, 73.63.Kv, 


INTRODUCTION 

At cryogenic temperatures semiconductor based quan¬ 
tum dots (QD) can be used to create well defined quan¬ 
tum states of arbitrarily few localized electrons. The elec¬ 
tron spins of these states provide a playground for explor¬ 
ing quantum mechanics in an interacting solid state envi¬ 
ronment and are heavily studied for possible applications 
in quantum information processing mm The coherent 
dynamics of electron spins can be accessed in an electron 
spin resonance (ESR) experiment. To control a QD based 
spin qubit on a time scale shorter than its dephasing time 
such an ESR experiment would require a magnetic held 
modulated at radio frequencies (rf) with an amplitude of 
few millitesla. Combining such a large rf modulation to 
an (externally applied) macroscopic magnetic held with 
cryogenic temperatures of T IK, required for long 
spin lifetimes, is a major technical challenge. Obstacles 
are oscillating strong mechanical forces between macro¬ 
scopic perpendicular magnets and eddy currents caused 
by induction, both causing severe heating and mechan¬ 
ical oscillations. To overcome these problems, on-chip 
methods to locally manipulate electron spins have been 
developed. A hrst breakthrough in locally controlling 
QD based spin qubits was based on the exchange cou¬ 
pling between two electrons located in adjacent tunnel 
coupled QDs [2l. This all-electrical method makes use 
of the direct dependence of the singlet-triplet splitting 
on gate voltages, while the latter can be rf modulated 
in a straightforward way mrnmm. Because the ex¬ 
change coupling between two electrons is subject to fluc¬ 
tuations of the local potential, it is, however, desirable to 


be able to manipulate the spin of a single electron local¬ 
ized in a QD as well. In a standard ESR approach this is, 
in principle, possible with an on-chip magnetic antenna 
M. This approach has the disadvantage of needing a 
relatively strong current through an on-chip micro-wire 
which causes parasitic heating of the sample. The ca¬ 
pacitive coupling between the strongly driven antenna 
and the QD leads can furthermore cause electron pump¬ 
ing via an unwanted modulation of the leads chemical 
potentials. Alternative methods are based on electric- 
dipole-induced spin resonance (EDSR) where a periodic 
spatial motion of an electron gives rise to an oscillat¬ 
ing (effective) magnetic field, an rf driving force. The 
rf spatial oscillation of an electron is thereby induced by 
modulating the voltage of one of the metal gates defining 
the QD. It has been demonstrated that the necessary in¬ 
homogeneous effective magnetic field can be provided by 
the spin-orbit interaction |H1 nu or even the spatial fluc¬ 
tuations of the hyperfine interaction between the electron 
and many nuclei [7 . Unfortunately, both these interac¬ 
tions also promote dephasing of the qubit; the spin-orbit 
interaction via coupling electrons and phonons, while the 
hyperfine interaction couples the electron spin dynam¬ 
ics to the thermal fluctuations of nuclear spins EDI m]. 
Consequently, it would be beneficial for spin qubit ap¬ 
plications to use materials combining a small spin-orbit 
interaction with no nuclear spins, e.g. 28 Si and 12 C. How¬ 
ever, this would require another mechanism to facilitate 
EDSR. 

An elegant option employs the inhomogeneous stray 
field near the edge of an on-chip magnet. A spatial os¬ 
cillation of an electron localized in such an inhomoge- 
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FIG. 1: Sample layout: (a) Scanning electron microscope im¬ 
age of the wafer surface. The GaAs surface is dark gray, gold 
gates are shown in yellow, magnetic cobalt gates in blue. Both 
magnets are ~ 2 pm long and ~ 60 nm high, the left one (L) 
is ~ 100 nm and the right one (R) ~ 230 nm wide. Red filled 
circles indicate possible QD positions in the two-dimensional 
electron system 85 nm beneath the surface; actual positions 
depend on gate voltages and disorder potential. The volt¬ 
ages applied to gates R and ~ are radio frequency modulated 
for EDSR measurements, (b) Magnetic force microscopy mea¬ 
surement of the magnets before the gold gates were processed; 
height profile in the upper panel and out of plane magnetiza¬ 
tion in the lower panel indicating single domain magnetization 
of both magnets. 


neous field then directly translates into a modulation of 
the magnetic field. In past experiments, relatively wide 
(width of ~ pm) on-chip magnets in the vicinity of dou¬ 
ble QDs have been used in order to create a strong field 
gradient 0 h uni na na m nu. The disadvantage of 
such a large magnet are its multiple magnetic domains at 
zero external magnetic field, B ext , which lead to a small 
and rather uncontrolled stray field. A sizable B ext on the 
order of a Tesla is then needed to align the domains and 
thereby create a strong inhomogeneous magnetic field at 
the QD. Multiple domains can be avoided by reducing 
the on-chip magnet’s lateral dimensions until its shape 
anisotropy yields a single-domain groundstate. In a pre¬ 
vious project, we have already realized an on-chip single¬ 
domain nanomagnet. It yields a sizable inhomogeneous 
stray field B nm independent of B ext and, therefore, pro¬ 
vides interesting possibilities for nanoelectronic circuits, 
in particular at B ext ~ 0. As an example we have demon¬ 
strated that this new regime can be utilized for very ef¬ 
ficient hyperfine induced nuclear spin manipulation and 
have indeed reached much stronger nuclear spin polariza¬ 
tions than previously reported for lateral QDs [22] . 

Here we present an innovative double QD hybrid design 
which incorporates two single domain nanomagnets. We 
replaced two of the usual gold gates with ferromagnetic 
cobalt gates (Fig. |TJl) . At small B ext , the two magnets 
can be magnetized in a parallel (as in Fig. Ef) or antipar¬ 
allel configuration, giving rise to two very different in¬ 
homogeneous magnetic field distributions, an interesting 
possibility for spintronics applications. The double QD 
is defined in the two-dimensional electron system (sheet 
density: 1.19x 10 11 cm -2 , mobility: 0.36x 10 6 cm 2 /Vs) of 
a GaAs/AlGaAs heterostructure 85 nm beneath its sur¬ 


face (Fig. So- We have prepared the double QD in the 
two-electron Pauli-spin blockade regime with one electron 
in each dot, in order to employ spin-to-charge conver¬ 
sion. To determine static properties such as the coercive 
fields of the two magnets we have used dc measurements 
and have explored the electron spin dynamics with EDSR 
measurements. 

Depending on the double QD geometry we have found 
either one or two electron-spin resonances. Two reso¬ 
nances corresponding to different B nm in the two dots 
would allow to study coupled spin qubits in a double QD. 
However, two resonances can only be resolved under three 
conditions: (i) a sizable magnetic field difference between 
the dots, (ii) a sufficiently large magnetic field gradient 
in each dot and (iii) a strong enough capacitive coupling 
between each dot and an rf-driven gate. While (iii) is 
straightforward to fulfill, in this article we demonstrate 
that the remaining conditions (i) and (ii) can be met 
by employing two single domain nanomagnets. As we 
merely replace gold gates by magnetic cobalt gates our 
scenario can be scaled up to multi qubit systems. Because 
single domain magnets are also useful at B ext = 0 they 
allow spintronics experiments beyond the scope of pre¬ 
vious experiments with only one (usually multidomain) 
on-chip magnet. 


EXPERIMENTAL SETUP 

Our measurements probe the dc current I (green ar¬ 
row in Fig. [lji) which passes through the double QD in 
response to a constant voltage V = {fi r — = 1 mV 

applied across it. Figure [2^i illustrates the double QD 
configuration by sketching the chemical potentials /iL,R 
of the leads and those of the relevant double QD states 
as horizontal lines, while vertical lines indicate tunnel 
barriers. We denote (n,m) the charge configuration of 
the double QD with n electrons in the left and m in 
the right dot and consider the single electron charge 
transfer characterized by the following tunneling cycle: 
(1,0) (1,1) -A (2,0) (1,0), where the transition 

(1,1) —)> (2,0) constitutes a bottleneck: Both configu¬ 
rations, (1,1) and (2,0), are composed of three triplets, 
collectively denoted by Tn and T 20 , and one singlet, Sn 
and S 20 • We define the detuning e as the energy dif¬ 
ference between Sn and S 2 o- For e < 0 tunneling pro¬ 
cesses (1,1) —)> (2,0) are blocked by energy conservation 
(Coulomb blockade) and / is close to zero. The exchange 
splitting between singlets and triplets is much higher if 
two electrons are in the same dot, J 20 Jin Hence, 
at e ~ 0 the T 2 o states are highly elevated compared to 
Tn. As a consequence, for 0 < e < J 20 the transition 
T 20 —>• Sn is forbidden by the Pauli principle (Pauli- 
spin blockade), as long as Sn and Tn remain decou¬ 
pled. In our case the inhomogeneous B nm mixes Sn and 
Tn states near where their eigenenergies are equal [22] . 
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FIG. 2: (a) Sketch showing the relevant two-electron levels of the double QD in Pauli-spin blockade and the lead chemical 

potentials (/xr — /xl = ImeV). All horizontal lines depict chemical potentials, vertical lines indicate tunnel barriers, blue areas 
correspond to occupied states of the degenerate Fermi liquid in the leads, (b) Corresponding charge stability diagram at 
B ext = 0, presenting the current through the double QD. The numbers of electrons charging the left, right dot in stable areas 
of the stability diagram are indicated in parentheses. For our EDSR measurements the system was pulsed between the two red 
dots; details in main text, (c) Current at e ~ 100 peV as function of B ex t which was swept at a rate of 30 mT/min. Arrows 
indicate the sweep directions, (d) Numerically calculated eigenenergies of the (1,1) states of the system Hamiltonian ([!]) as 
a function of B ex t- I^m) are almost identical to T± and 1^2, 3 ) are superpositions of To and S 11 . For the used parameters, 
e = 100 peV and t c — 4 peV, the (2,0) state has a much lower energy. Vertical steps at B ex t — 50, 80 mT indicate discontinuities 
of Bnm caused by the numerically included switching of the nanomagnets, (e) Overlap of the (1,1) eigenstates shown in d with 
the S 20 state (which would be zero for the Tr states in a homogeneous field): all four states in the inset. The main panel is 
a magnified view of the gray area. It contains the two states with the smallest overlap, the current limiting bottleneck states 
(almost T±). (f) Sketch of the four current resonances belonging to the four different magnets’ configurations. The red and 
blue line resemble the switching behaviour of the actual current in (c) (arrows indicate sweep directions). 


In the stability diagram plotted in Fig. IT this coupling 
gives rise to a narrow stripe of \I\ > 0 near e = 0. We 
stress that the hyperfine interaction which also couples 
Sn and T 11 , is only a weak perturbation compared to 
the effect of B nm . For e > J 20 the transition Tr —>> T 20 
lifts the Pauli-spin blockade and a sizable current flows 
(in our double QD J 20 — 300 peV). The tiny but non¬ 
vanishing current visible in Fig. It for 0 < e < J 20 is 
dominantly caused by higher order processes such as the 
co-tunneling Tr —>> T 20 —» (1,0), where T 20 is energeti¬ 
cally forbidden. 

In the stability diagram in Fig. [2J), the current carrying 
region corresponding to the absence of Coulomb block¬ 
ade is composed of two overlapping triangles. Above we 
have described the tunneling cycle which gives rise to 
the lower left triangle; the second triangle corresponds 
to an alternative cycle (2,1) -T (1,1) -T (2,0) -A (2,1). 
Nevertheless, the above explanations apply to both cy¬ 
cles, as they are both bottlenecked by the same transition 
(1,1)->(2,0). 

In previous devices on-chip magnets were separated 
from the heterostructure by a layer of metal gates and, 


with the exception of Ref. 22. in addition by a second 
electrically isolating layer. Here, we simplify the struc¬ 
ture and bring the magnets closer to the QDs by replacing 
two gold-gates (yellow in Fig.[]Ji), used to define the dou¬ 
ble QD, by ferromagnetic cobalt-gates (blue). Based on 
simulations with OOMMF [23] we have tailored the stray 
fields of the magnets and have optimized their geometries 
and positions to maximize stray field and field gradient 
between the two dots on the one hand and to guarantee 
full tunability of the double QD by applying gate voltages 
on the other hand. The advantage of using two instead of 
just one magnet is twofold: first, two magnets can be po¬ 
sitioned to provide strongly inhomogeneous and different 
magnetic fields in two adjacent dots (conditions (i) and 
(ii) above) which facilitates EDSR measurements in both 
dots. Second, at moderate B ext two separate magnets al¬ 
low for two very different stray field distributions across 
the double QD corresponding to either parallel or anti¬ 
parallel magnetization of the two magnets (see Fig. [ 5 ]). 
The magnetization of each nanomagnet can thereby be 
reversed by sweeping B ext beyond its coercive field and 
antiparallel to its present magnetization. The different 
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width of the two magnets gives rise to individual coer¬ 
cive fields. Consequently, we can choose between parallel 
and anti-parallel magnetization at relatively small B ext . 

SYSTEM HAMILTONIAN 

To model the dynamics of our double QD we assume 
that an electron localized in the left respective right dot 
experiences the local magnetic field B L,R = B ext + B R ^. 
We thereby neglect the hyperfine interaction between the 
electron and nuclear spins, the spin-orbit interaction and 


the exchange interaction which in our case are all small 
perturbations compared to the coupling induced by the 
inhomogeneous B nm . For simplicity we define the av¬ 
erage field in the two dots B = (B * l + B R ) /2, their 
difference field AB = B L — B R and the field operators 
B± = B x db iB y , A B± = (A B x ± iAB y ) /2 akin to spin 
raising and lowering operators. With the quantization 
axis z defined parallel to B ext , the matrix representa¬ 
tion of the (semiclassical) total Hamiltonian in the basis 
spanned by the diabatic singlet and triplet states {T + , 
T 0 , T_, Sn, S 2 o} is then 


T+ 

To 

T_ 

Sn 

S20 


( B z 

B_/V2 

0 

-AB_/-v/2 

0 ) 


t+ = in) 

B+/V2 

0 

B-/V 2 

AB z /2 

0 


To = (|n> + |lt»/V2 

0 

B+/V2 

-Bz 

A B+/V2 

0 


T_ = |U) 

-A B+/V2 

AB z /2 

AB-/y/2 

0 

tc/Zgm 


S 11 = (|tl)-|lt))/v / 2 

\ 0 

0 

0 

t c /2gii B 

-t/gn b ) 


s 20 = |o,tl> 


where t c denotes the interdot tunnel coupling between 
the two dots. The matrix representation 0 illustrates 
that the x- and //-components of the difference field, 
Amix T± with Sn , while the ^-component A B z 
mixes T 0 (which has no spin component along the z- 
axis). The average field B yields the Zeeman splitting 
of the spin-up versus spin-down states. Note that the 
off-diagonal terms B±, which mix T± with To, vanish if 
the quantization axis is chosen parallel to B instead of 

B ex t • 

Hyperfine and spin-orbit interaction would both con¬ 
tribute to various matrix elements including the singlet- 
triplet coupling constants. In our case, however, the lat¬ 
ter are far dominated by the time independent differ¬ 
ence field (and we formally neglect the former contribu¬ 
tions). In this way the nanomagnets provide a stabiliza¬ 
tion mechanism for appropriate qubit implementations 
which could increase the qubit coherence time in spite of 
the presence of nuclear spins or spin-orbit interaction. To 
fully determine our double QD hybrid system we need to 
know the nanomagnet’s strayfield as a function of B ext 
at the position of the two dots, the interdot tunnel cou¬ 
pling t c and the Lande g-factor inside the dots g. In the 
following we will employ dc current measurements and 
EDSR experiments to achieve this goal. 


DIRECT CURRENT MEASUREMENTS 

To experimentally determine the coercive fields of 
our nanomagnets we have measured the leakage current 

I through the spin-blockaded double QD while slowly 
sweeping B ext at constant detuning, e ~ 100 peV. The 
current at this configuration is sensitive to the mixing 
of the singlet and triplet states [3] and can be used to 
detect changes of the magnetic field differences between 
the QDs (22 : . 

In such a sweep experiment the current /(B ex t) might 
be influenced by dynamic nuclear spin polarization 
(DNSP) which can give rise to hysteresis as a function 
of the sweep direction of B ext [22j [24] [25]. However, it is 
also possible to avoid DNSP effects by preparing a fixed 
point at very weak polarization [22]. Here and in the 
EDSR experiments discussed below, DNSP is negligible 
and the apparent hysteresis of the measured /(B ext ) vis¬ 
ible in Fig. 0 has a different reason: it is related to the 
four distinct configurations of the magnets, each of which 
can be magnetized parallel or anti-parallel to B ext . 

In Fig. [ 2 ]^ we present /(£> ext ) for two sweeps in op¬ 
posite directions (dB/dB eXb = ±30mT/min). We have 
started the sweeps at large |T ext | to ensure that both 
magnets are magnetized parallel to B ex t- The current 
maxima near B ext = 0 occur where all three Tn states 
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are close to resonance with Sn. This is a consequence 
of Pauli-spin blockade where the leakage current is gov¬ 
erned by the singlet-triplet couplings: the Tn triplets 
mix with Sn and because Sn is tunnel coupled to the 
other singlet, S 2 o, the Tn triplets mix also with S 2 o- 
These mixings are strongest near the mutual resonances 
between Sn and Tn and zero far away from the corre¬ 
sponding resonances. As |B ext | is increased the T + and 
T_ triplets are more and more detuned from Sn, their 
mixing with S 20 also decreases, their decay (Tn S 20 ) 
slows down and they become the current limiting bot¬ 
tleneck states. Consequently, the current decreases. To 
illustrate this connection we numerically diagonalized the 
Hamiltonian in Eq. 0 and plot in Fig.[2]l the energies of 
the four relevant (1,1) eigenstates versus H ext for a sweep 
from negative to positive fields. The apparent avoided 
crossing at L> ext — 20 mT marks the point of minimal 
B ex t + B nm , where B z = 0. This field coincides with the 
current maximum for dB/dB ext > 0 (blue in Fig. [ 2 J 3 ) , 
because here the T±-S 2 o mixing has its maximum. For 
dB/dB e ^t < 0 the nanomagnets would be magnetized in 
the opposite direction and the current maximum would 
occur at L> ext ~ — 20 mT (as observed in the according 
measurement, red in Fig. [ 2 ]:). 

Continuing the sweep, we further increase B ext beyond 
the respective coercive fields of the two magnets where 
their magnetizations reverse and become again parallel 
to B ext . This change of magnetization instantly rear¬ 
ranges the overall magnetic field at the QDs and causes 
a steplike characteristic of the eigenenergies at the co¬ 
ercive fields. We show below that because of the direct 
relation between the singlet-triplet mixing and the nano- 
magnets’ configuration, this leads to the sudden changes 
of the measured current observed in Fig. [2J:, where we 
find coercive fields at H ext ~ 50 mT for the wider and 
-B ex t — 80 mT for the narrower magnet. These coercive 
fields are included in the numerics of Fig. [2]i . 

We remark that the observed current jumps occur very 
abruptly as a function of B ext . This underlines that the 
nano magnets are single domain and the single domain 
switches as a whole once the coercive field is reached. 

To phenomenologically explain the current /(L> ext ) in 
Fig. [2]^, we assume four distinct current maxima corre¬ 
sponding to the four possible configurations of our nano¬ 
magnets. By sweeping B ext we can switch between these 
configurations which causes the actual current to jump 
between the four maxima at the corresponding coercive 
fields. This simple, yet reasonable model is displayed in 
Fig. 0 , where we plot two pairs of Lorentzians (gray) 
reflecting the symmetry properties of the problem. The 
larger maxima correspond to the parallel and the smaller 
ones to the anti-parallel magnets’ configurations. For 
better comparability with our measurements we have 
added two colored curves which mimic the actual cur¬ 
rent jumps between the Lorentzians. The identical curves 
are also shown in Fig. 0 as black lines, where they 


reveal good agreement with the measured data. The 
two Lorentzians describing the parallel or anti-parallel 
magnets’ configurations, respectively, are equal in am¬ 
plitude and width. Most interesting is the observation 
that the overall current through the QD is considerably 
smaller if the two magnets are polarized anti-parallel to 
each other compared to their parallel configurations. It 
suggests that the anti-parallel magnetization causes a 
smaller singlet-triplet mixing of the bottleneck triplets 
T± than the parallel configuration. From our Hamil¬ 
tonian in Eq. 0 we see that the coupling between the 
T±-states and the singlet sub-space is proportional to the 
difference of the magnetic field component perpendicular 
to the quantization axis (approximately parallel to B ex t) 
between the two QDs. Thus it suggests that the per¬ 
pendicular component of the field difference between the 
two QDs is smaller when the magnetization of the mag¬ 
nets is anti-parallel compared to the parallel configura¬ 
tions. To check this, we have approximately determined 
the location of the two QDs (as depicted in Fig. [ 5 ] be¬ 
low) taking into account the numerically calculated B nm 
combined with results of the EDSR measurements (dis¬ 
cussed below), the gate voltage configuration (referred to 
as configuration I) and the measurement in Fig. [ 2 J 3 . This 
information provides the magnetic field components in 
our Hamiltonian in Eq. 0 (see table [T| and allows us 
to calculate its eigenenergies, shown in Fig. [2]l, as well 
as the mixings between the four ( 1 , 1 ) states and 820 - 
The latter are visualized in the inset of Fig. 0 as func¬ 
tion of L> ext swept from negative towards positive fields. 
The main panel is a magnification showing the mixings 
of the two bottleneck states T± near the current steps 
in Fig. 0. It also shows a steplike characteristic at the 
coercive fields of the two magnets and is considerably re¬ 
duced for the magnets in their anti-parallel configuration, 
namely between the two coercive fields. This strengthens 
the Pauli-spin blockade and explains the reduced current 
for the anti-parallel configuration of the nanomagnets in 
between the current steps) in Fig. [ 2 J 3 . We close this sec¬ 
tion by noting that an accurate prediction of the current 
which depends on B nm , B ext , e and t c would require a de¬ 
tailed density matrix calculation which goes beyond the 
scope of this article. 


ELECTRIC-DIPOLE-INDUCED SPIN 
RESONANCE MEASUREMENTS 

Mixing between any two singlet and triplet states is 
strongly enhanced where their eigenenergies are nearly 
equal. If detuned from this resonance, it is possible to 
actively drive transitions between two levels in an EDSR 
experiment which regains the resonance condition by ap¬ 
plying a proper rf magnetic field. We have performed 
our EDSR measurements in the (1,1) configuration with 
e ~ —100t c (lower red dot in Fig. [2 ]d) where the two elec- 
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FIG. 3: (a) ED SR measurement scheme, from left to right: 
initialization in Pauli-spin blockade —> separation of the QD 
states from the leads (Coulomb blockade) —>> EDSR manipu¬ 
lation by rf-modulation of gate voltage —> read-out (/ > 0 for 
EDSR resonance) and re-initialization, (b) Current I(B ext , f) 
through the DQD (bottom panel) and position of current 
maxima at EDSR resonance (top panel) while driving with 
the pulse sequence shown in (a), (c) Same as (b) but high res¬ 
olution measurement near the coercive field of the left magnet 
B z ^ 80 mT, where the resonance line forms a step. 


trons are strongly localized in the two respective dots. 
Consequently, we expect to find two distinct EDSR reso¬ 
nances at the respective Zeeman energies in the two dots 

hf = |<7MbB l,r | ~ g(i B (Bext + , (2) 

where / is the modulation frequency, h the Planck con¬ 
stant and the z-component of |Bj^|. The ap¬ 

proximation in Eq. © is fair for (B ext + B^\ z ) 2 > 
B ext and B r ^| 2 can have identical 
or opposite signs depending on the magnets configura¬ 
tion. The resonance condition Eq. © allows us to di¬ 
rectly probe the ^-factor as well as the z-component of 
B nm in the two dots and therefore also A B z . 

Our experimental EDSR sequence is sketched in 
Fig. [3ji. We start at e > 0 in the Pauli-spin blockade (up¬ 
per red dot in Fig. J2|)), which initializes the double QD 
with equal probabilities in one of the bottleneck states 
T + or T_. (The other two (1,1) states, T 0 and Sn, de¬ 
cay quickly if occupied and eventually the system stalls 
in T + or T_.) After 800ns we isolate the two electrons 
by pulsing the double QD deep into Coulomb blockade 
to e ^ — 100t c (lower red dot in Fig. 2b) by changing the 
gate voltages V& and W (see Fig. [lr) within ~ 2 ns. To 
avoid pulse transients effects we next wait 200 ns before 
we modulate W for 500 ns with a sine wave which causes 
both electrons to oscillate in real space. Due to the in¬ 
homogeneous B nm this rf modulation directly translates 
into oscillations of both B L and B R . Finally, we pulse 


back to our starting point at e > 0 for read out. If the 
rf modulation in both dots is off-resonant the double QD 
stays in Pauli-spin blockade and no current flows. How¬ 
ever, if the resonance condition Eq. © is fulfilled for one 
of the two electrons during the rf modulation, the Sn 
singlet state becomes occupied with a finite probability. 
As a consequence the Pauli-spin blockade is lifted dur¬ 
ing read-out. Performing a steady state measurement 
by periodically repeating this sequence at a frequency of 
~ 670 kHz we then measure a small leakage current. 

Typical results of such measurements are presented in 
Figs. ©b, c), where the leakage current is plotted as func¬ 
tion of H ext and /. At T> ext ~ 0 the Pauli-spin block¬ 
ade is lifted even without applying an rf-modulation as 
already seen in Fig. © and discussed there. This ef¬ 
fect gives rise to the broad frequency independent cur¬ 
rent maximum at H ext ~ 0. The rf modulation, how¬ 
ever, generates additional sharp but weak current max¬ 
ima along straight lines, where the resonance condition 
in Eq. © is fulfilled. Whenever B ext causes one of the 
nanomagnets to reverse its magnetization the resonance 
frequency suddenly increases according to the increase of 
the Zeeman energy. In Fig. ©3, T>ext was stepped from 
negative towards positive fields and hence the nanomag¬ 
nets reverse their magnetizations at the positive fields, 
H ex t — 50, 80 mT. Figure © shows a high resolution 
measurement of part of Fig. ©• It reveals one of the 
expected jumps in resonance frequency at H ext ~ 80 mT. 
Linear fits according to Eq. © suggest \g\ = 0.36 ± 0.01 
and B^ m \ z = 19 ± 1 mT if both magnets are magnetized 
parallel to B ext and B^ m \ z = 12 db 1 mT if the narrower 
magnet (with the larger coercive field) is aligned anti¬ 
parallel to B ext . 

We close the discussion of Figs. [3jb-c) with two re¬ 
marks: first, the switching of the large magnet cannot 
be observed in this measurement because it is masked 
by the broad current maximum around H ext = 0. Sec¬ 
ond, in Figs.[3jb-c) we find only one EDSR resonance for 
this particular gate voltage configuration (configuration 
I). We will show below that it belongs to the left QD and 
under which conditions a second EDSR resonance can be 
observed. 

In Fig. [4] we present a second EDSR measurement af¬ 
ter retuning the double QD such that each QD is situ¬ 
ated in close proximity to one of the two magnets. In 
this gate voltage configuration (II) the anti-parallel mag¬ 
netization of the nanomagnets did not result in clear 
EDSR resonances because of strong effects of DNSP at 
the required weak B ext values. Fortunately, the paral¬ 
lel magnetization of the nanomagnets allows for EDSR 
measurements at higher magnetic fields where DNSP is 
weak. In this regime we find two distinct EDSR reso¬ 
nances as can be seen in the raw data in Fig. [4^i plotting 
/(H ext ,/), in panel (b) indicating the positions of sharp 
current maxima in (a) and also in panel (c) presenting a 
single frequency trace at constant B ext . The latter con- 
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Driving frequency (GHz) Drivin 9 fre 9 uenc y f < GHz > 


FIG. 4: Coexistence of two EDSR resonances: (a) /(T ex t,/) 
through the DQD as in Figs. |3jb,c) but for a different gate 
voltage configuration (configuration II). (b) Position of the 
current maxima at EDSR resonances in (a). Lines indicate 
fits with eq. ©• (c) Typical current trace at constant B ex t 
along dashed line in (a). 



z (nm) 


z (nm) 


FIG. 5: OOMMF simulations of B nm [ 23 ] • Magnetic field 
component B n m|z along T ex t (z-axis) for (a) parallel and (b) 
anti-parallel magnetization of the two nanomagnets. B nm |z 
is relevant for the EDSR resonance condition, Eq. ©• (c, d) 
Absolute value of the gradient of the perpendicular compo¬ 
nent Bn m within the plane of the 2DES; B nm = B nm | z + Bn m - 
The EDSR signal strength scales with | V x , z Bn m | • Colored cir¬ 
cles in (a-d): Approximate QD center positions for the two 
different gate voltage tunings. Arrows indicate nanomagnets 
polarization directions. 


tains two sharp current maxima clearly indicating two 
distinct EDSR resonances. Fitting Eq. © yields the 
same \g\ = 0.36 =b 0.02 as above, but \ z = 9 db 2mT 
and B^ m \ z = 55 =b 3mT. The assignment of the EDSR 
resonances to either the left (L) or the right (R) QD is 
thereby based on our OOMMF simulations of B nm . In 
Figs.[5|a,b) we plot the B nm | 2 component which is rele¬ 
vant for the T± Sn transitions considered here for the 
two cases, a parallel and an anti-parallel magnetization 
of the nanomagnets. 

Comparing our numerical results with our EDSR mea¬ 
surements we find the approximate center coordinates of 
the QDs, which are marked as red circles in Figs. [5ja- 
d). The used gate voltages are in agreement with these 
locations, where details also depend on the disorder po¬ 
tential. The yellow circles in Figs. [5|a-d) indicate the 
approximate positions of the QDs for the previous gate 


voltage configuration (I) discussed above. 

The strength of each EDSR resonance (i.e. the height 
of each current maximum in Fig. should scale with 
the absolute value of the local gradient of the magnetic 
field component perpendicular to B ext , i.e. B^ m , along 
each electron oscillation path caused by the rf modula¬ 
tion. As we do not know the exact pathways we plot 
in Figs.[5jc,d) the absolute value of the two-dimensional 
derivative within the plane of the 2DES, |V x ^B^- m |. The 
red circles in Figs. [5jc,d) are clearly near derivative ex¬ 
trema while the situation is not so clear for the yellow 
circles. This observation provides a possible explanation 
for the missing second EDSR resonance in our first gate 
voltage configuration (yellow circles). There we can as¬ 
sign the observed resonance to the left QD as the sim¬ 
ulated B nm at the position of the left QD fits to the 
EDSR results for both magnet configurations but that 
of the right QD does not. Our numerical results at the 
marked QD positions are summarized and compared to 
our experimental findings in Table [l| Interestingly, for the 
parallel magnetization of the nanomagnets the difference 
fields A B x and A B y are quite weak compared to AB Z , 
while A B x and A B y are also sizable for the anti-parallel 
magnetization. This implies that the dynamics of the To 
state is quite different for the two cases, as the coupling 
between To and Sn is mediated by A B x and A B y . For 
instance, anti-parallel magnetization would be the better 
choice for defining a qubit based on the states To and Sn, 
while the parallel magnetization would be a good choice 
for a qubit based on T± and Sn. 


configuration I (single EDSR resonance) - T nm (mT) 


magnet iz. 

QD 

Bnm | x 

B nm | y 

Bnm \z 

A B x 

A By 

A B z 

=4 

L O 

RO 

-15 

-12 

6 

-1 

17 (19) 

22 

-3 

7 

5 


L O 

RO 

i i 

CO to 

D- kO 

i i 

9 (12) 

17 

1 

2 
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configuration II (two EDSR resonances) - B nm (mT) 


magnet iz. 

QD 

Bnm \ x 

T nm | y 

Bnm | z 

A B x 

A By 

A B z 

=4 

L • 

R# 

14 

13 

19 

26 

10 (9) 
50 (55) 

-1 

7 

40 (46) 

* 

L • 

R# 

-10 

8 

-11 

29 

10 

48 

18 

30 

38 


TABLE I: Magnetic field components generated by the two 
nanomagnets (for parallel [=4] vs. anti-parallel magneti¬ 
zation) at the approximate QD positions marked in Fig. [5] 
by corresponding circles. Field values are calculated with 
OOMMF [23] ; field strengths derived from measured EDSR 
resonances in parentheses. 





















































External magnetic field B ext (mT) External magnetic field B ext (mT) 


FIG. 6: Double quantum dot with one instead of two single 
domain nanomagnets reproduced from Ref. m- (a) Scan¬ 
ning electron microscope image; color coding analogue to EK 
(b) OOMMF simulations of the perpendicular field gradient 
along B ext aligned with the z-axis. This is also the main di¬ 
rection of the rf electron displacement. The lower panel is a 
horizontal cross section of the color plot at x — 100 nm. (c) 
/(Bext, /) through the DQD as in Figs. showing a sin¬ 

gle EDSR resonance at a modulation amplitude of ~ 15 mV. 
(d) Multiple EDSR resonances at about twice the modulation 
power. 


SAMPLES WITH ONE NANO MAGNET 

In the following, we discuss an EDSR measurement of 
a sample of our previous generation of double QD designs 
[26] . It contains just one single domain nanomagnet (see 
Fig. [6^i) located on top of a gold gate. The voltage on 
the same gate is modulated to drive EDSR. This design 
is especially simple as the magnet axis, the rf-modulated 
gate and the symmetry axis of the QD (z-axis) coincide. 
It is justified to assume that the gate voltage modulation 
entails a motion of the QD electrons mostly also along the 
z-axis. Simulations of the nanomagnet shown in Fig. [6)3 
reveal a sizable magnetic field gradient at the left QD 
while it almost vanishes at the right QD at a larger dis¬ 
tance to the magnet. Our measurements agree with the 
simulations and reveal an average field of approximately 
|B nm | ~ 22 mT within the left QD and a coercive field 
of the nanomagnet of 52 ± 2 mT. As in our EDSR ex¬ 
periments, the rf-magnetic field modulation is produced 
by driving an electron along a time-independent slanting 
magnetic field, only electrons in the left QD can be ma¬ 
nipulated by EDSR. Consequently, in this sample only 
a single EDSR resonance corresponding to the left QD 
is expected. This is in contrast to our sample with two 
nanomagnets, which can be tuned such that a sizable 
field gradient exists in two QDs giving rise to two EDSR 
resonances associated with two separate QDs (see Fig.[4|. 
As expected, we find only a single EDSR line in the ex¬ 


periment in Fig. [6]3. The g-factor is identical to the one of 
our previous sample, g = 0.36 ±0.01, while both samples 
feature similar QDs based on the same wafer. 

In Fig. [6]i we present a similar EDSR measurement 
as in Fig. [6 }j, but with about twice the modulation 
power. We observe a transition from a single resonance 
at huo = g/iB^ext small modulation powers in Fig. & 
to multiple resonances at nftw = g/ip^ext; n = 1,2 ,... 
for larger modulation powers in Fig. Such a be¬ 

havior can be explained in terms of higher order har¬ 
monics generation and has also been observed to such a 
high order in comparable EDSR experiments in an InAs 
nanowire based double QD m and to second order in 
another GaAs based double QD [28] . A related example 
of higher order harmonics generation are Landau-Zener- 
Stiickelberg-Majorana interference oscillations ES ng. 
The multiple resonances here (with variable slopes) are 
fundamentally different from the two resonances (with 
identical slopes) observed in our previous sample with 
two nanomagnets. There, they are caused by the differ¬ 
ent magnetic fields in two QDs and, in contrast to the 
higher harmonics, observable for low driving powers. 

CONCLUSIONS 

In summary, we have explored a hybrid nanostructure 
consisting of a double QD incorporating two single do¬ 
main nanomagnets with different coercive fields. The 
magnets’ properties agree well with numerical simula¬ 
tions. By sweeping an external field it is possible to 
reverse the magnet polarizations one-by-one and to di¬ 
rectly measure their coercive fields. The magnetization 
of each nanomagnet switches almost instantly at its co¬ 
ercive field as expected for a single domain magnet at 
very low temperature. Each switching event modifies the 
local magnetic field distribution and gives rise to a dis¬ 
tinct current jump in a dc transport measurement in the 
Pauli-spin blockade regime. In a radio frequency EDSR 
experiment the switching of a magnet generates a shift 
of the resonance frequency. Compared to the larger mul¬ 
tiple domain magnets our single domain magnets gener¬ 
ate sizable field gradients even at zero external field and 
therefore allow experiments in a regime where the relative 
field difference between adjacent QDs is stronger than the 
average field. In contrast to magnets with multiple do¬ 
mains, the single domain nature guarantees a stable field 
distribution over a large range of external field values. 
The disadvantage of a somewhat smaller field gradient 
due to the smaller size of the magnet can be compen¬ 
sated by using multiple single domain nanomagnets. The 
combination of several magnets provides a higher degree 
of control of the overall field distribution as the coercive 
fields of each magnet can be predetermined by design. In 
summary, coupled QDs including multiple single domain 
nanomagnets represent a promising approach for future 
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spin qubit circuits desired for quantum information or 
related spintronics applications. 
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APPENDIX 


Methods: Magnetic field simulations were carried out 
with the 3D solver oxsii 1.2a5 of the OOMMF toolkit 
[2E[ . A simulation gridsize of 5nm (in plane) and 2.5 nm 
(out of plane) was chosen, comparable to the magnetic 
exchange length of ^ 3.5 nm in cobalt. The exchange 
stiffness of cobalt, 30 x 10 -12 J/m, as well as the satura¬ 
tion magnetization, 1.4 x 10 6 A/m, enter the simulation 
as external parameters. We included a reasonable 10° 
correction between the 2DES-plane and the direction of 
B ext which provides the best fit to the measured sample 
characteristics. 
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